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Abstract 

We perform lattice calculations of the spatial 't Hooft fc-string tensions, 5^, in the decon- 
fined phase of SU(iV) gauge theories for N = 2,3, 4, 6. These equal (up to a factor of T) the 
surface tensions of the domain walls between the corresponding (Euclidean) deconfined phases. 
For T > T c our results match on to the known perturbative result, which exhibits Casimir 
Scaling, &k oc k(N — k). At lower T the coupling becomes stronger and, not surprisingly, our 
calculations show large deviations from the perturbative T-dependence. Despite this we find 
that the behaviour oc k(N — k) persists very accurately down to temperatures very close 
to T c . Thus the Casimir Scaling of the 't Hooft tension appears to be a 'universal' feature 
that is more general than its appearance in the low order high-T perturbative calculation. We 
observe the 'wetting' of these fc-walls at T ~ T c and the (almost inevitable) 'perfect wetting' 
of the k = N/2 domain wall. Our calculations show that as T — > T c the magnitude of 5&(T) 
decreases rapidly. This suggests the existence of a (would-be) 't Hooft string condensation 
transition at some temperature which is close to but below T c . We speculate on the 'dual' 
relationship between this and the (would-be) confining string condensation at the Hagedorn 
temperature Tjj that is close to but above T c . 



1 Introduction 



In the Euclidean formulation of SU(iV) gauge theories at finite temperature T, deconfinement 
at T = T c is associated with the spontaneous breaking of a centre symmetry. This implies 
the existence of iV degenerate phases for T > T c in our Euclidean box. When such phases 
co-exist they are separated by domain walls. The tension of these domain walls is equal, up to 
a factor of T, to the spatial 't Hooft string tension. This is an interesting dynamical quantity 
that, for T < T c , is in principle related to confinement. 

As T — > oo one has g 2 (T) — > jl] and so one can calculate the spatial 't Hooft tension in 
perturbation theory P] El- These calculations [2] show that this tension satisfies Casimir 
Scaling (to 2 loops). This is an intriguing result given that the usual confining string tension 
between static sources is also known to satisfy Casimir Scaling to quite a good approximation 
[HI El IE] ■ Clearly it would be interesting to know whether the high T Casimir Scaling that is 
characteristic of low order perturbation theory, persists to lower values of T where the coupling 
becomes strong and non-perturbative effects should be important. In particular it would be 
interesting to know what happens at T ~ T c where one can imagine making contact with the 
confining phase. 

In this paper we answer this question using numerical lattice techniques applied to SU(2), 
SU(3), SU(4) and SU(6) gauge theories. The plan of the paper is as follows. In Section I2~TI 
we discuss in more detail the theoretical background to this problem. We then provide a 
summary of our lattice setup in Section 12.21 In Section 12.31 we describe how we can impose 
the existence of a domain wall by using twisted boundary conditions and then, in Section |23J 
we describe how we can calculate (the derivative of) the surface tension from the average 
action. In Section 12.51 we give the high-T prediction for the profile of the domain wall and 
describe how we obtain this profile in our lattice calculation. We complete Section El with 
a discussion of finite volume corrections in Section 12.61 and of 'wetting' in Section 12.71 We 
then move on in Section El to our calculations. Section 13.11 contains our results on the /c-wall 
tension and makes the comparison with perturbation theory. In Section I3~2l we calculate the 
domain wall profiles and again compare with perturbation theory. In Section 13.31 we look for 
the appearance of a layer of confining phase in the centre of the /c-wall as T — > T c , pointing 
out that one can use the pattern of 'wetting' to tell us something about the domain wall 
between confining and deconfining phases. We finish in Section 13.41 with a conjecture about 
the existence of a dual to the usual string condensation 'Hagedorn' transition. Section HI gives 
our conclusions. 

While this work was in progress a paper appeared that addresses some of the same 
questions that we do, using a novel lattice technique that allows the direct calculation of the 
tension (rather than its derivative, as in this paper). The calculations, for SU(4) at a fixed 
value of the lattice spacing, observe consistency with Casimir Scaling down to T = 1.2T C just 
as we do. However our calculations test the Casimir scaling formula more tightly, because we 
work with SU(6) as well as SU(4), and we vary the lattice spacing so that we can address the 
continuum limit to some extent. In addition our range of T extends much closer to T c . 
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2 Preliminaries 



2.1 Background 

As we remarked above, deconfinement in SU(iV) gauge theories is associated with the sponta- 
neous breaking of a Zn centre symmetry. The symmetry arises from the fact that the partition 
function is unchanged if we impose that the time torus is periodic only up to a global gauge 
transformation belonging to the centre of the group, z k = exp{2iirk/N} G Z N . (The group 
element z k contains a factor of the N x N unit matrix which we have suppressed. In fact it 
will be convenient to use the same notation, z k , both for the phase factor and for the corre- 
sponding matrix. Which of the two we mean should always be clear from the context.) The 
obvious order parameter is the trace of the Polyakov loop, l p , which is a closed Wilson line 
that winds once around the Euclidean time torus and whose trace therefore acquires a factor 
of z k under this symmetry transformation. In the confined phase the effective potential for 
the Polyakov loop trace averaged over the volume, l p , has its minimum at l p — and so the 
symmetry is explicit. In the deconfined phase the minimum is at l p ^0 and so the symmetry 
tells us that there are N degenerate minima: the symmetry is spontaneously broken and there 
are N possible deconfined phases. If l p oc z k we label the deconfined phase by k. When two 
of these deconfined phases, k\ and k 2 say, co-exist they will be separated by a domain wall 
whose surface tension will depend on k = ki — k 2 , as well as on N and T. At high T the 
coupling on the relevant length scale, g 2 (T), becomes small and one can calculate the domain 
wall surface tension in perturbation theory. To two loops one finds [TJ E] : 

k(N - k) ^7Sm {] - '^ (T)N) (1) 

where c 2 — 0.09. (While the precise value of c 2 depends on the coupling scheme used [T], there 
appears to be only a ~ 10% variation between physically reasonable schemes.) The recently 
calculated 3-loop contribution [2] turns out to be negligible and so we do not include it here. 
The spatial 't Hooft string tension (see below) is related to the wall tension by 

5* = (2) 

(We use this notation because a k usually denotes the Wilson-loop /c-string tension.) 

The factor of T 3 in eqn(£Q) follows from dimensional analysis once T is large enough for it 
to be the only relevant scale, and to two loops a^r is only a function of the 't Hooft coupling 
A = g 2 N [HE]. The leading depend ence on the coupling is 1/g rather than 1/g 2 because there 
are no semiclassical domain walls: they are stabilised at one-loop. Although one might worry 
about the infrared divergences that eventually obstruct the perturbation expansion, numerical 
tests in 2+1 dimensions J3, have shown that the lowest order perturbative calculations are 
remarkably accurate at all temperatures as long as they are not very close to T c . Of course, 
as we reduce T we increase g 2 (T) and at some point the perturbative expansion becomes 
unreliable. For example if we use the 1-loop perturbative expression for the Debye screening 
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mass, m D = T g 2 (T)N/3 [4 , and equate it to the lattice values calculated in [5] then we 
would find 

g 2 {T)N -7.5 , T ~ 2T C (3) 

at which point the 2-loop correction in eqnQ becomes ~ 60% of the 1-loop value and it is 
clear that we are past the point where we can be confident in perturbation theory (and indeed 
in the reliability of the calculation of m^). At these values of T numerical calculations are 
the only unambiguous means of determining the physics. 

The factor k(N — k) in eqn(Q) is the k dependence of the Casimir, Trj{T a T a , when the 
represention 1Z is the totally antisymmetric representation of a product of k fundamentals of 
SU(iV). Thus it is the factor one obtains when calculating the Coulomb interaction between 
sources in such a representation. In .0=1+1 SU(iV) gauge theories, the tension of the confining 
fc-string that emanates from such a source has precisely this dependence. In fact there are 
old speculations [S] that this 'Casimir Scaling' holds in D=3+l and numerical calculations 
show that it is in fact a good (but not exact) approximation to such fc-string tensions [7 j 
- and an even better one in D=2+l A question we will try to answer in this paper is 
whether the Casimir Scaling in eqn(^), which is an outcome of the low-order perturbative 
calculation, survives at very much lower T where perturbation theory becomes a manifestly 
poor approximation. This non-perturbative question will be addressed using numerical lattice 
techniques. 

What connects this question to the role of Casimir Scaling in confinement, is that these 
domain walls are closely related to the centre vortices that provide a possible mechanism for 
confinement fH\ \1'2\. Such a vortex will disorder any Wilson loop that encircles it, and the 
same for appropriate pairs of Polyakov loops. (We temporarily locate our discussion in the 
in the more easily visualisable case of D=2+l [Sj.) At low T, in the confined phase, one 
postulates a condensate of such vortices, and a simple statistical calculation (see e.g. [13]) 
shows how this leads to linear confinement. Since the Wilson loops are in a Euclidean space- 
time plane, the vortex cross-section will be in the same plane. As the Euclidean time extent is 
reduced so as to increase T, at some point the vortex core is squeezed by the short temporal 
direction and becomes periodic in time. At this point one can no longer find a Wilson loop that 
encircle the vortex, but it will still disorder Polyakov loops that it separates. This co-incides 
with the onset of the deconfined phase, which is a 'Higgs' phase for these vortices, in which 
they become real, 'massive' objects. They are now precisely the domain walls which separate 
deconfined phases characterised by different Polyakov loop expectation values. The domain 
walls/vortices we discuss wind around a spatial torus. However they could equally well be of 
finite length, ending on appropriate monopoles [Hj. Since the vortex cross-section is in 
a space-time plane, these monopoles will have world-lines in an orthogonal spatial direction 
(returning to -D=3+l). They form the boundary of a spatial 't Hooft loop, and the domain 
wall tension is just the corresponding string tension (multiplied by T). Where this tension is 
non-zero, e.g. for T > T c , these loops are not condensed. If however their tensions exhibit 
Casimir Scaling all the way down to T c , where they might condense, then this may provide 
an extra ingredient in our (hopefully) growing understanding of confinement. 

As an aside, we remark that although these domain walls certainly feature in the 4- 
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dimensional Euclidean field theory, they almost certainly do not entail real domain walls 
in the gauge theory at finite T. This has no implications for the calculations described herein. 
(A physical manifestation of these walls is described at the end of Section l2~3*l ) 



2.2 Lattice setup 

We discretise Euclidean space-time to a periodic hypercubic lattice. The lattice spacing is a 
and the lattice size, in lattice units, is L M in the \i = 0, 1, 2, 3 directions. We choose the \i = 
direction as the short Euclidean time direction that generates our temperature: 

T — — j— • (4) 

We will typically choose L\ = L 2 and place the domain wall to span the L\ x L 2 torus. 
Typically L3 3> L\ since it will need to contain 1 or 2 domain walls separating regions in two 
phases. In addition, for thermodynamics to be applicable, we need aLi,aL 2 3> 1/T. 

The fields are SU(iV) matrices, Ui, on each link I of the lattice. (We sometimes label Ui 
more explicitly as U^{n) where [i is the direction of the link and n = {no, n} is the site from 
which it emanates.) We label the elementary squares (plaquettes) of the lattice by p and the 
ordered product of link matrices around the boundary of p by U p . The partition function is 
then 

Z=[HdU l e^ s (5) 
J 1 

and we use the plaquette action 

S = E^ 1 - ^ ReTrf U ( 6 ) 

p 

Analysing the continuum limit of this path integral, one finds that 

o 2N 2N 

9 9l( a ) 

where we have used the fact that g 2 is a lattice bare coupling to write it as some lattice scheme 
running coupling, g^, on the length scale a. To approach the N —>■ 00 limit smoothly, we need 
to keep fixed the 't Hooft coupling 

= g\l)N (8) 

(where I is fixed in units of some well behaved quantity such as the mass gap) and so it is 
useful to define the 't Hooft inverse coupling 

7 = ^- = ^— (9) 
<-2N? gl(a)N- W 

If we keep 7 fixed as we take N —>■ 00 then we ensure that a remains fixed in physical units 
for large enough N. (See Fig.l of |3] for an explicit lattice demonstration of this.) 
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2.3 Imposing domain walls 

Our order parameter will be based on the Polyakov loop 

l p (n) = -Trl[u (n ,n). (10) 

"0 = 1 

Consider the transformation where one multiplies all timelike link matrices at some fixed 
n = n' Q by a non-trivial element of the centre of SU(N), z k = e 2m N G Z N : 

U (n = n' ,n) — > z k U (n = n' ,n) Wn. (11) 

A plaquette involving one such link matrix will necessarily involve the conjugate of a second 
and the centre factors will cancel, z k z k ^ = 1 so the plaquette remains unchanged. Since the 
measure is also unchanged, such a field has the same probability as the original field and we 
have a corresponding Zn symmetry. A winding operator such as the Polyakov loop in eqn fjlOj) 
will change under the symmetry transformation: l p — > z k l p . It therefore provides an order 
parameter for the spontaneous breaking of the symmetry. In the confined phase where (l p ) = 
the symmetry is explicit, but in the deconfmed phase where (l p )^0 it is spontaneously broken 
with N degenerate phases in which (l p ) = z k c(f3) where c(/3) is a real-valued renormalisation 
factor. Such phases can co-exist at any T > T c and in that case they will be separated by 
domain walls. If the order parameter in the two phases differs by a factor z k we shall refer 
to the domain wall as a A;- wall. The domain wall may, for example, extend right across the 
1-torus and 2-torus, and its world volume will wrap around the time torus. In such a case 
periodicity in £3 means there must be at least two walls. Since the Polyakov loops on either 
side of the domain wall differ by some factor z k e Zn, we see that the wall disorders Polyakov 
loops. Although we cannot encircle the wall with a Wilson loop in the (3,0)-plane, since the 
wall wraps right around the time torus, if we could then the Wilson loop would acquire a 
factor z k as well. Thus, as remarked earlier, the domain wall is nothing but a spatial 't Hooft 
disordering loop in the (l,2)-plane, that has been squashed in the short temporal direction. 
We describe it as squashed because the size in the x 3 -direction, where it is not squashed, is 
of the order of the inverse Debye mass 

i Ts °°^S»f^° <i2) 

to leading order in perturbation theory 0j. 

To study a fc-wall we 'twist' the action in eqnfJUJ) so as to enforce the presence of a single 
domain wall. The twisted action is defined by 

^^^(l-^ReTr^p)^}). (13) 

p 

where z(p) = 1 for all plaquettes except 

z{p = {[iv, Xf,}) = z k = e 2m ^ (iv = 03; x = x ', x 3 = x 3 ', x h2 = 1, L lf2 . (14) 
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That is to say, the plaquettes in the (0,3)-plane that are at the fixed values x = x ' and 
x 3 = X3 and at all values of Xi and x 2 are multiplied by z k G Z N . This is a plane of plaquettes 
that wraps around the 1 and 2-tori at xq = xq and £3 = X3. Although the plane appears to 
be in a particular location this is an illusion. We can move it and deform it by redefinitions 
U\ — > ZkUi of suitable subsets of links. It is clear from considering such redefinitions that the 
Polyakov loops on either side of such a plane of twisted plaquettes will differ by a factor of z^. 
Periodicity in x 3 then demands that at some £3 the Polyakov loop must suffer a compensating 
factor of zjJ so as to return to its original value. This requires an interface - a domain wall. 
Thus using such a twisted action we ensure that in the deconfined phase each configuration 
possesses at least one domain wall. This provides a convenient way to study different fc-walls 
corresponding to different choices of k in eqn (jl3|) . 

We will not discuss the interpretation of these walls in any detail, and instead refer the 
reader to [SJ for a more explicit discussion in the readily visualisable context of D—2+1, 
and to ^1 for a more general theoretical review. Here we content ourselves with a few remarks 
for orientation. A finite line of twisted plaquettes enforces the presence of a Dirac string so 
that we have monopoles at the ends of the line The return flux between the monopoles 
is a centre vortex. A plane of such plaquettes thus corresponds to a 'dual' Wilson loop where 
a centre monopole propagates around the perimeter of the loop. There is now a vortex sheet 
that also ends on the perimeter of the loop. In our construction the loop is spacelike and we 
then extend the plane of twisted plaquettes right around the corresponding spatial two-torus, 
so that the boundary disappears and the vortex sheet, which now wraps around the two-torus, 
becomes decoupled in position from the centre Dirac sheet. Nonetheless the presence of the 
latter enforces the presence of the former. One may regard the plane of twisted plaquettes as 
effectively constituting an ultraviolet /c-wall at zero cost. In the picture where vortex sheets 
are degrees of freedom relevant to confinement in the Euclidean theory at T = ^B], the 
idea is that a closed fc-vortex sheet can 'thread' a Wilson loop (just as a vortex loop would 
do in D = 2 + 1) and the Wilson loop then acquires a factor z k - A condensate of such loops 
then leads to an area decay and linear confinement. These vortex sheets will in general close 
upon themselves rather than around the (infinite) torus. Our periodic construction is a device 
to impose the presence of a vortex sheet so that we can study it. In the deconfined phase 
such a spacelike vortex sheet is squashed by the short Euclidean time direction so that it also 
extends around the time torus. And simultaneously the vortex sheet becomes real just as a 't 
Hooft-Polyakov monopole becomes real in the Higgs phase of the Georgi-Glashow model. It 
is in this phase that we can study it in the hope that what we learn will eventually have some 
bearing on confinement. 

Finally we remark that this whole discussion becomes more transparent if we perform a 
cyclic relabelling of our co-ordinates, x Q — > x\ etc., so that we have a sytem at low T = 
1/0X3 <C T c but on a spatial three-torus where one of the tori is very short. In this case our 
Zn monopole sources are timelike and the dual 't Hooft loop is in a space-time plane. The 't 
Hooft loop is imposed by the sheet of twisted plaquettes which generate the monopole Dirac 
sheet. It is now the short spatial torus that squashes the return flux between the monopoles. 
If we make the sheet of twisted plaquettes periodic in both directions then the squashed return 
flux manifests itself as a domain wall separating two phases in which the Polyakov loop around 
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the short spatial torus has non-zero vacuum expectation values that differ by a factor of Zf. 
between the two phases. This relabelled description has the advantage that it corresponds to 
real T ~ physics on such a spatial volume, in contrast to the original description which (for 
subtle reasons) does not imply the existence of real domain walls in high-T gauge theories. 



2.4 Domain wall tension 

Let us calculate the average action with and without a fc-twist as defined above. Then, on a 
lattice that is of fixed size in lattice units, 

AS k = (S k ) - (S ) - — — - - (15) 

where is the surface tension (energy per unit area) of the domain wall and A = a 2 L\L2 
is its area. We assume the spatial dimensions are large enough that finite volume corrections 
are negligible. (We will of course need to check that this is so for our particular calculations.) 
We see from eqnflTJ that we expect 

AC 8c k w A T ^ ,. MN-k) fc' 2 L,L 2 8 1 

ASk = ap— = a(u) ^rrwi-wT0m (16) 

when T = l/a(/3)L is large enough that leading order perturbation theory is accurate. The 
additional factor of a(Lo) contains the 0(a 2 T 2 ) = 0(1/ Lq) lattice correction to the continuum 
formula in eqn(^). We have obtained it by a numerical evaluation of the (non-closed) analytic 
expression given in [3] and the values are listed in Tabled for a range of aT. The calculations 
in this paper will be for L = 4, 5 where the correction is modest. 

To compare the perturbative expression in eqn (flfi|) with our calculated values of AS k 
we need to know the coupling g 2 (T) in terms of g\(a) and hence f3. The simplest choice is 
obtained by asserting that to leading order the scale and distinction between running couplings 

is irrelevant, so we can use g 2 (T) = g\(a) and then 

AS t t-oo ,. . 2jt 2 k(N-k) 1 1 

ur 2 = a(Lo, vi^^ii7F (17) 

In practice this is too crude. It is well known that g\(a) is a bad running coupling in that 
it has very large higher order corrections, so that it is only at unattainably asymptotic T 
that we could be confident of eqn(JT7J) holding. The simplest, if partial, remedy is mean field 
improvement ^7j 

■N " gfta 
If we replace g 2 (T) in eqn (J16)) by gj(a) then we obtain 

2n 2 k(N — k) 1 1 

- ^WtNllWt (19) 



Pj = (3 x (-Tr[/ P ) = (18) 
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where the difference between the second line and the first is 0(1/(3*) and is therefore negligible 
even to 2-loop order. (If we use the expansion to 0([3~ s ) obtained in JH] for the average 
plaquette, we see that even at our maximum value of a, corresponding to a = 1/4T C , the 
numerical correction is no more than ~ 17%.) In a similar way we obtain from eqn(|T|) 

AS k 2 ^° P A^x{l + 4l^i- (20) 

We do not have a calculation of the lattice correction for the 2-loop term so we use no correction 
in that case. The error from this should be small. The expressions in eqns (jl9l20j) are what 
we shall usually use when making comparisons between perturbation theory and our lattice 
results. 

There still remains the ambiguity associated with the scale T at which gj should be eval- 
uated. Although, as we remarked above, the scale is irrelevant at one loop, this is really a 
formal statement. The point is that with a reasonable coupling, using an appropriate scale 
is the best way to ensure that higher order corrections are likely to be small. Of course, for 
L = 4, 5 plausible scales such as irT and 2ttT are numerically close to 1/a, so it may be that 
using gj(a) for g 2 (T) is reasonable. However, it would be useful to have some estimate of 
how eqn(|19|) changes if we use the scale aL , say, rather than a in g 2 . Near T c this is a large 
spatial scale and so we need to use a lattice coupling scheme that we can have real confidence 
in. For this purpose we shall use the Schrodinger functional running coupling defined and 
calculated in [TH]. We shall leave the detailed discussion and calculation to the appropriate 
point in Section El 



2.5 Domain wall profile 

In addition to calculating the surface tension of the /c-wall as described above, we can also 
calculate how the average Polyakov loop varies as one passes from one vacuum to the other 
through the domain wall, just as was done for the D=2+l SU(2) case in 

Suppose the /c-wall spans the (1,2) torus and is centered at z = where z £ (— oo, +oo) 
and we use continuum variables. Suppose the Polyakov loop satisfies 

l p (z = -oo) = l ; l p (z = +oo) = e 2 ™^. (21) 

Then to leading order in perturbation theory one finds |1 3 

l p (z) = ^Tre 2 ™^ (22) 

where is a diagonal matrix with {Yk]u = k for i = 1, ...,N — k and {Yfc}™ = k — N for 
% = N — k + 1, ...,N and q(z) is a function that is independent of k. It is obtained from a 
one-loop effective potential calculation for the Polyakov loop PjJ Ej which leads to 



e V4^ 

-Tz 
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Eqns (j22l23j) provide the 1-loop perturbative result for the domain wall profile for all k and N. 

To compare this with our lattice calculations we have to choose a coupling g 2 (T) and, just 
as for the surface tension, we choose to use the coupling gj(a) defined in eqn ([T8|) . In addition 
the average Polyakov loop in the deconfined phase does not have unit length, but rather 

\(l p )\=c(P). (24) 

This can be interpreted as the self-energy correction to the static source and is dominated by 
ultra-violet fluctuations which are insensitive to which phase we are in. Thus if we want to 
compare to the numerically calculated Polyakov loop we can simply renormalise the predicted 
Polyakov loop by a constant factor at all z 

l p (z) -> c{(5%{z) (25) 

where c(/3) can be determined from field configurations that are in a pure phase (and will 
depend on iV as well as on 0). In principle we could obtain lattice corrections to eqns (j22|23j) 
analogous to the factor a(Lo) for the surface tension. However we expect these to be small 
and so we neglect them. Thus eqns ()22J23pi8J2oj) provide the 1-loop perturbative result for the 
domain wall profile against which we shall compare our numerical lattice results. 

2.6 Finite volume corrections 

In principle the thermodynamic limit demands that we take Lj — > oo ; % = 1, 2, 3 at any fixed 
L and (3. In practice what this means is that we need to make these lengths large enough so 
that finite V corrections are small. We start with the most obvious. 

2.6.1 wall fluctuations 

The fluctuations of the domain wall lead to corrections that are geometric and will not possess 
the k, N dependence of the domain wall tension - in exactly the same way as the Luscher 
correction to the mass of a confining string does not depend on the representation of the 
sources. Thus we want to make the wall area, L\ x L 2 large enough for these corrections to 
be negligible. The leading corrections are typically |20] 

c 1 

L0L1L2 CT^y 

3x/3yW)iV Ll 
C 4tt 2 k(N-k) LiL 2 

0.01 (26) 

where we have used the fact that c ~ 1 (201, e Q n © as some guide to the largest value of 
g 2 (T)N and the fact that for our calculations we have Lq = 4, 5 and Li^ = 20. Thus with 
our choice of lattice sizes any such corrections will be very small. 
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2.6.2 wall depth 

We must also make sure that L3 is long enough to accommodate the full depth of the domain 
wall. As pointed out in P] this requires a much large L3 than one might naively expect. If we 
define the full width of the L3 = 00 wall, dw, as extending between the points at which the 
Polyakov loop reaches two-thirds of its vacuum value, then the domain wall starts expanding 
once aL 3 < 3dw and completely disappears once aL 3 ~ 2d w . We therefore explicitly monitor 
the domain wall profile in each calculation and make sure that L 3 is large enough that there 
are no significant finite-/^ corrections. 

2.6.3 spatial symmetry breaking 

This correction is the most subtle and is typically relevant to calculations at very high T. In 
our case the problem arose when we performed calculations in SU(4) on 4 x 12 2 x 40 and 
4 x 12 2 x 60 lattices at (3 = 18, 20 respectively. The latter (3 corresponds, very roughly, to 
T ~ 1000T C and we expected to obtain a good agreement with perturbation theory. However, 
while our result for the k = 2 domain wall was in agreement with eqn ()19)) . the k = 1 result 
was about twice too large and even larger than the k = 2 result, which was clearly unphysical. 

The cause of this phenomenon turns out to be a centre symmetry breaking in the spatial 
direction. That is to say, the average Polyakov loop in the /x = 1 or /x = 2 direction acquires a 
non-zero vacuum expectation value. In the case under discussion, the centre breaking occurred 
in the simulation with k = 1 twist but not for the untwisted k = case. Thus the difference 
in actions was anomalous. 

The practical cure, followed in the calculations in this paper, is obvious - calculate spatial 
Polyakov loops and choose the spatial lattice large enough to ensure that the symmetry 
breaking does not occur. For practical and theoretical reasons it would be interesting to 
understand this phenomenon sufficiently well to be able to estimate the spatial lattice size 
needed. This we now turn to. 

At high T, if we perform a Fourier decomposition in xq of the SU(iV) fields, we can keep 
just the lowest constant mode since there is an effective mass gap to the higher modes of 
0(tcT). The Euclidean time integral in the action just gives 1/T and the action becomes 
that of a D = 2 + 1 gauge theory coupled to an adjoint scalar field (the remnant of the Aq 
field). The coefficient in the action of the plaquette term is the corresponding D=2+l inverse 
coupling, which we call gf , and which clearly satisfies 



at tree level. (An approximation that can be systematically improved upon). If we now 
neglect the massive adjoint scalar field (again an approximation that can be improved upon) 
we obtain a D = 2 + 1 SU(iV) gauge theory on a L\ x L 2 x L 3 lattice with the usual plaquette 
action. All this is the well-known phenomenon of dimensional reduction at high T. (See [21] 
for a recent discussion.) 

Suppose for a moment that L 2 ,-L 3 ^> L\. Then as we increase (3% this D = 2 + 1 SU(iV) 
gauge theory will at some point deconfine and at this point the average Polyakov loop in 




(27) 
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the shorter x% direction (which plays the role of the Euclidean time direction in the D = 
2 + 1 system) acquires a non-zero expectation value. This is precisely the 'spatial symmetry 
breaking' in the parent D = 3 + 1 gauge theory that we discussed above. When will it occur? 
This is easy to estimate. Using eqn ([2Tj) we have a specific value of (3% corresponding to given 
values of (3 and Lq. The string tension, a 2 a, can be obtained for any N and for any [3$ from 
the calculations in [22] • (It helps to know the average plaquette for the interpolation.) We 
now note that in D = 2 + 1 one finds T c ~ yfa [23]. (This is based on SU(2) and SU(3); 
calculations for all N are in progress [211 •) Thus in the case cited above, with (3 = 20 and 
Lq = 4 and N = 4, we have (3$ ~ 80 which turns out to correspond to a^fo ~ 0.08 |2Sj so 
that the critical value of L\ will be Lf" rf = l/aT c ~ Xja^fa ~ 12. Since in D = 2 + 1 one has 
a^fa (x l//?3 at large fa, the value of Lf lt corresponding to (3 — 18 in the D=3+l theory will 
not be very different. Finally, the fact that L 3 is finite and that L 2 ^Li but rather L 2 = L\ 
means that the transition is smeared in f3 3 (and hence in f3) and that it may occur for one of 
the two transverse directions chosen at random. Thus the observed scenario where on some 
lattices the transition occurs and on others it does not, is just what one expects. Using the 
above approach one can estimate the critical value of L 12 for any particular (3 and L <C L\p 
and thus avoid the corresponding spatial symmetry breaking. 

It is interesting to note that at high T in D = 2 + 1 there should be no such spatial 
symmetry breaking because the dimensionally reduced D — 1 + 1 gauge theory confines at 
all T. (The confining string cannot vibrate in one space dimension.) Indeed careful finite V 
studies in D=2+l [3] have not observed any unexpected behaviour of the domain wall tension 
at very high T. It is perhaps over-familiarity with D = 2 + 1 that led us to be surprised by 
what occurs in D = 3 + 1. 



Since the phase transition is first order for N > 3 one can have the confined phase as well as 
the various deconfined phases co-existing at T = T c . We can always interpolate between any 
two deconfined phases by passing through the confined phase, and whether the system chooses 
to do so or not will depend on the relevant wall free energies, i.e surface tensions. A fc-wall 
between two deconfined phases will be stable against breaking up into two confined-deconfined 
domain walls, with surface tension a c d, provided that 



In this case it is only at some T < T c that the fc-wall would split. Another possibility is 



This corresponds to 'perfect wetting'. It has been conjectured [21] that this occurs in the 
SU(3) gauge theory, where there are only k — 1 walls. (See 0| for a recent discussion.) A 
third possibility is that the continuation from higher T would give 



2.7 Wetting at T ~ T c 



a k w {T c ) < 2a cd {T c ). 



(28) 



a k w {T c ) = 2a cd {T c ). 



(29) 



a k w (T c ) > 2a cd (T c ). 



(30) 
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In that case there would be some critical value of T above T c where a pair of confmed- 
deconfined domain walls would have the same free energy and from there onwards, as T — > T c , 
the fc-wall will be on the cr^(T) = 2a c d{T) branch. (So that in practice the condition in 
eqn(j2nj) will be satisfied.) 

We recall that the confined and deconfined phases correspond to the N + 1 degenerate 
minima of the effective action, calculated as a function of l p . The confined phase is at l p — 
and the N deconfined phases at l p = c((3)e 2mk l N : k = 0, N — 1. A domain wall between two 
phases is a tunnelling path between the corresponding minima that minimises the effective 
action. (A proper analysis, as in PQEj for very high T, needs to express the effective potential 
not just in terms of the trace of the Polyakov loop but in terms of its eigenvalues.) At high 
T it is enough to keep the effective potential and a term quadratic in derivatives QjJ |2] . In 
this case one has the usual simple intuition about how tunnelling paths are determined by 
the potential. As T = l/aL decreases the coupling g 2 (T) increases and eqn(J2HJ) tells us that 
the domain wall will become thinner, oc l/g(T), on the scale of 1/T. This could imply that 
higher order derivative terms will become more important in the effective action. In reality 
however, as we shall see, the walls become thicker rather than thinner as T — > T c , suggesting 
that the simple picture of a potential plus quadratic derivatives might remain reliable beyond 
its naive domain of validity. 

If we increase T slightly away from T c then the minimum at l p — rises and the confining 
phase becomes metastable. If we keep increasing T then eventually this minimum ceases to 
play any role and we need only consider the N deconfined minima. At high T one can calculate 
these tunnelling paths in perturbation theory PUI2] as summarised in Section l2*31 For example, 
one finds that the k = N/2 domain wall that runs between the k = and k = N/2 phases 
passes through J p = 0. In practice, as we shall see below, this behaviour persists all the way 
down to T = T c . Given this we expect a^ N ' 2 (T c ) = 2a c d(T c ); i.e 'perfect wetting'. Indeed 
we expect this equality to hold at any T close enough to T c for <r c d(T) to make sense. Thus 
we expect that at T = T c the k = N/2 wall will break up into two confined-deconfined walls. 
Since at T c the confined and deconfined free energies are equal, these two walls can separate 
and rejoin. 

As an aside we remark that the above scenario allows us to calculate T c on arbitrarily 
large volumes, in contrast to the standard calculations such as in j^J. For example, direct 
tunnelling at T = T c from the confined to the deconfined phase on a lattice larger than 
16 3 5 will simply not occur in a SU(6) Monte Carlo calculation of feasible length 0. The 
reason is the thermodynamic suppression of the large domain walls that occur as a necessary 
intermediate step in the tunnelling. Here, by contrast, our twisted boundary conditions enforce 
the presence of a fc-wall that can have enough free energy to split into two confined-deconfined 
walls irrespective of how large is the volume. In this context (3 C can be identified by the 
fact that the growth and shrinkage of the confined and deconfined sub-volumes follows an 
unbiased random walk. In principle this provides a practical technique for directly locating 
T c on arbitrarily large volumes. 

A ky^N/2 domain wall, on the other hand, does not pass through l p — at high T. Thus 
an alternative tunnelling path from the k = to K = N/2 minima is not via l p = 0, but rather 
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via the paths joining the sequence of minima K — — > 1 — > ... — > N/2. While this is clearly a 
local minimum of the effective action, this minimum turns out to be sufficiently higher than 
the direct path through l p = that it does not play a role in practice. More interesting is the 
reverse possibility: for tunnelling between the k and k' phases, there is a minimal path that 
first runs from the k minimum to l p = (just like the first half of the k = N/2 domain wall) 
and then runs out from l p = to the k' minimum. We would expect this to have an effective 
action equal (approximately) to that of the k = N/2 tunnelling path and thus to be too large 
to contribute at high T (where we know that we have paths that satisfy Casimir Scaling). 
However as T drops to and below T c there will certainly be a point at which the path through 
l p — comes to have a lower action than the path not through l v = 0. It could of course be 
that the shape of the effective potential also changes so that the latter path smoothly deforms 
into the former at T ~ T c . However our calculations (see below) suggest that what actually 
occurs is that there is a tunnelling from one kind of path to the other. So if Casimir Scaling 
persists all the way down to T c , we expect that it will only be at some T below T c that the 
k = 1 wall breaks up into two cd walls. In this sense we do not expect 'perfect wetting' for 
domain walls in SU(3), since these are necessarily k = 1. 

The above remarks suggest a possible way of determining whether a c d is O(N) or 0(N 2 ), 
something which direct calculations j3] have so far been unable to determine unambiguously. 
Since we have found that <7jy(T) satifies Casimir Scaling very accurately down to T ~ 1.02T C , 
it is plausible to assume that it will continue to do so down to T = T c . In that case the 
iV-dependence is <r^(T c ) oc k(N — k) and so 2o d = a^ N ^ 2 (T c ) cx N 2 . More generally by 
observing for which values of k and iV the domain walls wet at T = T c , one can imagine being 
able to pin down a cd quite accurately. The calculation in this paper does not have the range 
of k and iV to make this a realistic goal, but we can at least attempt to gauge the potential 
of such a strategy. 

3 Results 

The strategy of our calculation is as follows. We first perform calculations of k — 1 and k = 2 
domain walls in SU(4) at a very high temperature, T ~ 1000T C , where we can expect our 
results to match perturbation theory quite precisely. We then repeat the calculations much 
closer to T c , where perturbation theory is a priori no longer useful. We do so for L = 4 
and for L = 5, so that at each T we have results at two values of a. i.e. a = 1/4T and 
a = 1/5T. This allows us to make some statement about the continuum limit. We then 
perform calculations of k — 1, 2, 3 domain walls in SU(6), for L t = 5, which together with the 
SU(4) results provide us accurate information on the full k and iV dependence of the surface 
tension. As a check on the iV-dependence of the low T behaviour of the k = 1 surface tension 
we also perform some SU(2) and SU(3) calculations. 

We simultaneously calculate the domain wall profiles and compare them with perturbation 
theory. As T — > T c we look for the breakup of the /c-wall into two confined-deconfined domain 
walls separated by a region of confining phase ('wetting' and 'perfect wetting'). In principle 
this provides a way to obtain information about the iV-dependence of the corresponding 
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surface tension, a cd . 

Finally we list the values of /3 that correspond to T = T c for the calculations in this paper: 

{10.4865(3) - 0.067(7)/^T 3 : L = 4 ; 517(4) 
10.6373(6) - 0.090(17)/1/T 3 : L = 5 ; £17(4) (31) 
24.5139(24) - 0.112(19)/VT 3 : L = 5 ; SU(6) 

The last two values are taken from OES] while the first is a calculation performed in this paper, 
Exhibiting the volume dependence allows us to correct for a finite volume shift, although in 
practice this will be completely negligible for us. 

3.1 Surface tension 

The parameters of our high statistics calculations are listed in Table El and the corresponding 
values of the average plaquette are listed in Table 01 We have deliberately chosen the values 
of (3 so that the values of T/T c are closely matched across the various SU(iV) groups. (For T c 
we use the values in eqn ()31)l and for interpolations we use the the string tension values in [HI] 
as well as the extra calculations listed in Table EJ One can instead interpolate using the values 
of /3 c (Lq) in ESj and this will give values of T/T c that differ at the 1% level from the ones 
we list.) These values, T ~ 1.88T C and T ~ 1.02T C , take us very close to the deconfinement 
transition itself. In addition we have performed a calculation for SU(4) at a very large value 
of j3, which corresponds to T ~ 1000T C , where we expect to find good agreement with the 
perturbative result. 

In Table El we list the results of our calculations of ASifa. We also list the 2-loop pertur- 
bative expectation obtained from eqn()20|). using for g 2 (T) the mean- field improved coupling 
gj(a) defined in eqn (|T8*|) . 

Consider first the values of AS^ that we obtain in SU(4) at (3 = 20, which corresponds 
to T ~ 1000T C . We observe excellent agreement with perturbation theory for both k = 2 and 
k = 1 walls. The ratio of the k = 2 and k — 1 values therefore exhibits Casimir Scaling, i.e. 
is oc k(N — k). At such a very high T all this is to be expected - but it provides a reassuring 
check on our calculations. 

Moving on to the other SU(4) calculations with aT = 1/Lq = 1/4, we observe that as 
T —>T C the value of AS^ grows much more rapidly than the perturbative prediction. Indeed 
at T ~ 1.02T C the discrepancy is very large, a factor of ~ 18, and even at T ~ 1.88T C it is a 
substantial factor ~ 3. This is essentially telling us that da^/dT becomes much larger than 
the low-order perturbative expectation as T — > T c . Since at large T the domain wall tension 
tends to its perturbative value, what this implies is that cr^(T) is increasingly suppressed 
relative to its perturbative value, as we approach T c . Later on in this Section we shall provide 
an estimate of this suppression. 

Although we find that A Syy deviates strongly from low-order perturbation theory as T 
approaches T c , we also see from Table Elthat the ratio of AS^r 2 to AS^ 1 continues to satisfy 
Casimir Scaling, to a high precision, even down to T ~ 1.02T C . This is a remarkable and 
unexpected result. 
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All these calculations are at a fixed value of the lattice spacing in units of T, i.e. aT = 1/4. 
At this value of a we expect typical lattice spacing corrections to both the non-perturbative 
and perturbative calculations to be small but not negligible (see e.g. Table [l}. Nonetheless 
it would be useful to check that the effects we are seeing do not change significantly in the 
continuum limit. To obtain some control over that limit we repeat our calculations for Lq = 5, 
i.e. aT = 1/5. The leading lattice corrections should be 0(a 2 T 2 ) and so should be reduced by 
a factor ~ 2/3 when we go from Lq = 4 to Lq = 5. We see from Table El that at T ~ 1.88T C 
we observe a discrepancy with perturbation theory that is essentially identical to the one we 
saw at the coarser lattice spacing. Thus we can infer that this discrepancy is a continuum 
effect. For T ~ 1.02T C the ratio to the perturbative result is ~ 15 for Lq = 5 rather than 
the ~ 18 we saw for Lq = 4, but since this ratio is varying very rapidly with T when we are 
this close to T c , it is not clear whether this apparent variation with a is due to lattice spacing 
corrections or is due to a slight mismatch between the values of T/T c at the two values of a. 
In any case it is clear that a large and growing mismatch with perturbation theory as T — ► T c 
is a feature of the continuum theory. Nonetheless we see that the ratio of the AS^ continues 
to accurately satisfy Casimir Scaling when we decrease a, so that is also a property of the 
continuum theory. 

We now turn to our SU(6) calculations. These have been performed for Lq = 5 rather 
than Lq = 4 so as to avoid the first order bulk transition [5 j that separates the lattice strong 
and weak coupling regimes for N > 5. This and the fact that our statistics are smaller (the 
computational cost grows roughly oc iV 3 ) means that our results are significantly less precise 
than for SU(4). Nonetheless it is clear from Table El that we observe at both values of T 
precisely the same discrepancy with perturbation theory as we saw for SU(4) at the same 
value of Lq. In addition the AS*^ ratios continue to satisfy Casimir Scaling (albeit only at the 
2<T level at T ~ 1.02T C ). This and the fact that the discrepancy with perturbation theory is 
the same as for SU(4) tells us that the domain wall tension has no factors of k and N except 
for the Casimir scaling factor k(N — k) and its dependence on the 't Hooft coupling, g 2 N. 
(The latter follows from the fact that everything is only a function of T/T c and that fixed T c 
is obtained by keeping g 2 N fixed at large N [S].) 

Finally we note in Table that our calculations of the SU(2) and SU(3) k = 1 domain wall 
tensions show a very similar discrepancy with perturbation theory. (The SU(2) transition 
is second order and this prevents us from performing reliable domain wall calculations at 
T ~ 1.02T C .) Thus the suppression of cr^(T) as T — > T c is largely independent of N. 

To illustrate the above, we plot in Fig. [T] the values of AS^ 1 for our L = 4 SU(4) calcu- 
lations, and compare to the two-loop perturbative expectation, eqn(JT7j), using the mean-field 
improved coupling, as well as the unimproved lattice bare coupling (3. (The fact that both 
the lattice and perturbative calculations exhibit Casimir Scaling means that we would obtain 
a very similar comparison in SU(6).) We see the growing discrepancy as T — > T c that was 
discussed above. We also see a substantial difference between the two perturbative predic- 
tions. It is therefore relevant to ask how much of the observed discrepancy with perturbation 
theory might be due to the fact that our choice of coupling or scale is becoming inappropriate 
on the large distance scales associated with the limit T ~ T c . To address this question we 
shall choose some standard 'good' coupling whose running has been evaluated down to low 
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enough energy scales that we can extrapolate with some confidence to scales ~ T c . (This 
approach is at best qualitative: a proper calculation would require us to redo every part of 
the calculation in whatever coupling scheme we use, and this (more formidable) calculation 
we do not attempt.) An example of such a coupling is the Schrodinger functional coupling 
whose running has been calculated non-perturbatively in ^HJ- In Appendix A, we describe 
the evaluation of this coupling on the length scale T _1 = aL Q and we show the result in Fig|21 
Note that the direct calculation jTHj goes down to about T ~ 2T C and from there onwards 
we extrapolate using the '3-loop' beta function calculated in [TH]. We use this coupling to 
obtain AS^ and we show the result in Fig. [TJ We see that while using this coupling reduces 
the discrepancy at T ~ 1.88T C , it it is still far from large enough to bridge the observed gap 
at T ~ 1.02T C . This suggests that for T — > T c the domain wall tension becomes increasingly 
dominated by non-perturbative physics. 

Returning to the evidence for the Casimir Scaling of the fc-wall tensions, we average our 
accurate SU(4) ratios for aT = 0.20,0.25 giving 



A qk=2 



'1.327(26) 
1.345(27) 
1.341(49) 
1.333... 



T ~ 1.02T C 

T ~ 1.88T C , . 

32 

T ~ 1000T C 
Casimir Scaling. 



This (together with the results at T/T c = 2.3,1.5,1.2 of P) provides impressively accurate 
evidence for the Casimir scaling of domain wall tensions (or 't Hooft string tensions) from 
T = oo all the way down to T ~ 1.02T C . 

The fact that, when T decreases towards T c , da^/dT increases more strongly than ex- 
pected from perturbation theory, tells us that must decrease more rapidly than expected 
as T — > T c . In principle we can calculate function of T by interpolating <9a{y /d(3 in (3 

and then integrating from large (3, where it asymptotes to the perturbative value, down to the 
desired value of T. In practice our calculations are not dense enough in (3 to determine the 
interpolation unambiguously. (They were not intended for this purpose). Therefore the best 
we can do is to obtain a qualitative picture by assuming some functional form for da^/df] 
and fitting it to the few calculated values we have. Choosing the L = 4 SU(4) calculation, we 
display in Fig. [T] our best fit to ASw based on a simple modification of the one-loop formula 
to the fc-wall free energy: 

jpk 

n [1 loop] + aexp(-&Vft - Ac) (33) 



L ± L 2 T 



(The fit is also constrained to go through the very large (3 value which is not shown in the 
Figure.) Using this fit we then obtain the the T-dependence of a\y shown in Fig. El We 
display the uncertainty that arises from the errors in the fitted parameters, but it is clear that 
we cannot reliably estimate the systematic error inherent in the choice of fitting function. 
However we can see from Fig. ^ that our interpolation is conservative in that, if anything, it 
errs on the side of minimising the value of daw/dT in the region T > 1.02T C . We also show 
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Fig. El in the two- loop perturbative prediction (using the improved lattice coupling). The 
qualitative picture is that is strongly suppressed at T < 2T C , well below the perturbative 
expectations, to the point of being close to zero at T = T c . This is not inconsistent with our 
picture of these domain walls as being squashed but 'real' Zn vortices that condense into the 
vacuum in the confined phase. 

3.2 Profile 

As described above, imposing suitable twisted boundary conditions in the deconfined phase 
ensures that each lattice field contains a single domain wall. (Assuming, as is always the case 
for us, that the thermal production of extra pairs of domain walls is completely negligible.) 
If we average the Polyakov loop over the transverse co-ordinates we obtain the profile l p {x^). 
This profile maps the interval x% G [0, L — 1] to a line in the complex plane corresponding to 
the values that l p {x^) takes as x% runs over the X3-torus. This line can be directly compared to 
the perturbative one loop prediction in eqns (j21H23|) . Note that this comparison is independent 
of the choice of g 2 {T). (As we see from eqn(|23|) the value of the coupling determines the rate 
at which l p runs along the line in the complex plane, as £3 is varied, but not the actual locus 
of the line.) Thus any disagreement between calculated and perturbative loci cannot be put 
down to the use of an inappropriate coupling or scale. 

In Fig. Eland Fig. 03 we compare the one-loop perturbative prediction for the SU(4) k = 
1 locus, to what we obtain at T ~ 1.88T C and at T ~ 1.02T C respectively, in our lattice 
calculations. Plotting the fluctuating locus for each lattice field would clutter the plot so 
what we show is an average line. These results are also characteristic of what we observe for 
Lq = 5 and for SU(6) and other k. What we observe is not only excellent agreement at the 
very highest T (not shown here), but also, from Fig.EJ that the agreement remains remarkably 
good down to T ~ 1.88T C . This implies that the shape of the relevant region of the effective 
potential remains close to the one-loop perturbative form, even at these low T. By contrast, 
once we descend very close to T c , as in Fig. El at T ~ 1.02T C , the average locus has moved 
far from the one-loop curve suggesting that the effective action is now far from its high-T 
perturbative form. This observation makes it all the more remarkable that we continue to see 
Casimir Scaling at such very low T. 

In SU(4) the average k = 2 wall profile between the k = and k = 2 phases, must lie 
along the real axis by symmetry, just like the perturbative prediction. It is interesting to ask 
whether this is achieved by averaging over, say, two sets of tunnelling paths fluctuating around 
two degenerate minimum action paths, each away from the real axis, but which are related 
by Im/p — > — Im/p, or if there is only one such minimum tunnelling path, as in perturbation 
theory. In Fig. Inland in Fig. [7] we show the loci for individual lattice fields containing a k = 2 
wall, to illustrate that the latter is the case, both at T ~ 1.88T C and at T ~ 1.02T C . Thus 
we can expect, as discussed in Section 12.71 that the k = N/2 wall will indeed split into two 
confined-deconfined domain walls precisely at T = T c . 

It is also relevant to look at the corresponding plots for the k = 1 SU(4) wall. In Fig. [HJwe 
show what happens at T ~ 1.88T C . More interesting is T ~ 1.02T C . We have seen in Fig. El 
that at this very low T the average locus is deformed so as to pass not very far from l p = 0. 
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As we pointed out in Section [2~71 we expect that as T — ► T c there will be a second minimum 
action tunnelling path that passes along the real and imaginary axes near l p = 0, and which 
becomes the global minimum near T c . Does the average in Fig. El include some paths that are 
tracing this alternative trajectory? In fact this appears not to be the case, at T ~ 1.02T C , as 
we see from Fig. El In the next Section we shall see that at T = T c things do in fact become 
different. 

As we have just seen, in practice the fluctuations are small enough that we can identify 
quite accurately the profile and hence the centre of the domain wall in each lattice field, at any 
of the above values of T. We can then translate that centre to a common value of X3 and take 
an average over all the lattice fields so as to obtain the averaged domain wall profile. Let us 
relabel the position in the 3-direction by z, with the domain wall centered at z — (L— 1)/2 and 
z G [0, L — 1]. In Fig. El we plot the average of the real part of the Polyakov loop as a function 
of z for the L = 4 SU(4) calculation at T ~ 1.88T C . Note that for the k = 1 SU(4) wall 
the real part is not symmetric about the centre: one must exchange the real and imaginary 
parts. We do not show the imaginary part, to avoid cluttering the plot, but the same analysis 
and results follow in that case. We also show curves obtained using the one-loop formula in 
eqns f!2HI23|) with four different choices of coupling. Simply using the bare coupling g 2 = 2N/ j3 
does not reproduce the shape at all. The mean-field improved coupling defined in eqn ([TS|) 
works better but not very well. The Schrodinger functional coupling discussed in Appendix 
A works very well indeed. The last coupling, which is simply defined so as to minimise the 
discrepancy between the one-loop formula and the lattice result, works only slightly better 
than this. 

At T ~ 1.0 2T C we already know from Fig. that no coupling will work well. We illustrate 
this in Fig. ^2 where we plot the modulus of the Polyakov loop against z and compare it to 
the various perturbative predictions. Here we see the interesting feature that the 'better' the 
coupling, the worse is its description of the profile. This is in contrast to what we find in Fig.lT2l 
in a similar plot at T ~ 1.88T C . The reason for this is that the 'better' couplings are larger 
at lower T than the lattice bare coupling, and eqn (|23|) then predicts a thinner domain wall in 
lattice units (recalling that T = 1/aLo is fixed in lattice units). Here however the observed 
thickness of the wall actually grows as we lower T/T c from 1.88 to 1.02. This would appear to 
contradict asymptotic freedom which tells us that g 2 (T) grows as T j and hence, by eqn(j22J), 
that the wall narrows. This points to the fact that near T c low-order perturbation theory is 
qualitatively, and not just quantitatively, unreliable. It is interesting to note that precisely 
the same phenomenon is found [3] to occur for the lightest mass that couples (strongly) to 
the Polyakov loop. (At high T this is just twice the electric screening mass m^.) As we see 
from eqn (|12|) . we would expect mo to approach its high-T linear behaviour from above since 
g 2 (T) decreases with T; however lattice calculations show that the approach is from below 
Since I/mo is the natural length scale for the thickness of the domain wall, this observation 
ties in with what we have found here. 
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3.3 Wetting 



Motivated by the discussion in Section 1277] we have performed a set of calculations extremely 
close to T = T c . They are not intended to be accurate enough to provide a useful estimate of 
AS; rather we wish to see if the fe-wall separates into two cd walls at T ~ T c in a Monte Carlo 
sequence of several thousand field configurations. We perform such calculations for k = 1,2 
walls in SU(4) on a 4 x 100 x 20 2 lattice and for k = 1, 2, 3 walls in SU(6) on a 5 x 100 x 20 2 
lattice. We use long lattices so that if a fc-wall breaks up into two cd walls there is enough 
space for two well separated walls. 

Our observations are summarised in Table |B1 and Table The first thing we note is that 
the transition from no sign of wetting to the complete splitting of the domain walls takes place 
over a very small range of T which is precisely around T = T c . (See eqn([3ip.) For k^N/2 we 
first have a profile that is well away from the origin ('no'). Then at lower T we typically see 
fluctuations of the wall profile which deform it so that it sometimes passes through the origin 
('via origin') but does not visibly break up. If we lower T further the wall fluctuates more 
strongly, so that the profile not only goes through the origin, but the wall splits, reforms, and 
again reassumes a profile well away from the origin ('mixed'). Finally, at lower T still, the 
fc-wall splits completely ('splits'). All this takes place over an interval 8T/T C ~ 0.005. For 
the k = N/2 walls not all these stages are identifiable because the profile normally passes 
through the origin. But a precurser to the separation, a central 'flattening', can be identified 
and presumably overlaps with the 'via origin' and 'mixed' stages for k^N/2. 

It is interesting to ask what happens to the k = 1 SU(4) profile, shown in Fig. |H1 for 
T ~ 1.02T C , when we lower the temperature to T ~ T c . To illustrate this we show in Fig. [TBI a 
set of 10 profiles obtained over a sequence of 1000 Monte Carlo sweeps and in Fig. ITU over a 
later sequence of 1000 sweeps, all at T ~ 1.006T C . Each locus has been obtained by averaging 
over subsequence of 100 field configurations. We observe that the domain walls in Fig. 
have profiles that fluctuate around an average locus that is similar to the one we saw in Fig. |H| 
The profiles in Fig. HU on the other hand, fluctuate around a quite different locus that passes 
through l p = and which is identical to the minimum k = 2 path up to the 7r/2 rotation at 
the origin. What we appear to see is a tunnelling between these two paths and this strongly 
suggests that the effective potential has two minimal paths separated by a ridge through which 
the system has to tunnel. Once it has tunnelled to the set of paths passing through l p = it 
will readily split into two cd walls at T = T c . One would expect that the surface tension of 
the k = 1 and k = 2 walls that pass through l p = will be approximately the same rather 
than being related by Casimir Scaling. The statistics of our calculations near T = T c do not 
allow a direct test of this reasonable conjecture. 

We observe that, within our resolution, the wetting of the SU(4) k = 2 and k = 1 walls 
appears to occur simultaneously and it appears to occur at T = T c . This would appear to 
tell us that a^r(T c )^t2a c d(T c ) for both k = 1 and k = 2 at iV = 4. For SU(6) there is perhaps 
some evidence that the k = 2 and k = 3 walls begin to split slightly before the k = 1 wall. If 
this hint were to survive a more careful calculation, then we would have the beginnings of a 
calculation of the iV-dependence of a c d as oulined earlier. 

The fact is that all the k- walls appear to break up at almost the same value of T again 
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raises the question whether this can be consistent with the fc-wall tensions being different and 
related by Casimir Scaling as T — > T c , since in that case it is clear that one would expect the 
higher-/c walls to split at a higher T than the lower-fc walls. 

To address this question and also to understand why everything occurs in such a very 
narrow range of T around T c , we need to look at the magnitudes of the quantities involved. 
First consider the 'natural' magnitude of near T = T c as given by the one loop part of 
eqn((TJ) and using a realistic coupling. Plugging in the value that the Schrodinger functional 
coupling takes at T c in SU(3), g(T c ) ~ 2.2 and assuming g 2 (T c )N is constant, we obtain 

4-tt 2 T 3 

4,(1 - loop) = k(N - k)— ^___ ~ 2k(N - k)Tl (34) 

We now compare it to the latent heat [S] 

|^0.35AT 2 T c 3 + O(^). (35) 

It is natural to compare L h divided by T c to because the latter is a free energy per unit 
area that is effectively integrated over the thickness of the wall, which is 0(1/T C ). We observe 
from eqn ()34|) that for the maximal wall, k = N/2, cr{y(l — loop) ~ 0.5iV 2 T 2 and this is of the 
same order as the value of L^/T c . This seems reasonable, given that there is one scale in the 
problem. Now consider a C( i- Although, as we have emphasised, it is not known if it increases 
as N 2 , such a fit is certainly possible, and one then finds [3] 

a cd ~ (0.014iV 2 - 0.104)T C 3 . (36) 

This is very small compared to the previous two quantities; e.g. for SU(4) it is about a factor 
of 50 smaller than Lh/T c . The reason for this mismatch between these two energy scales is 
not entirely clear (but see below). Nonetheless, as we shall now argue, it is what leads to the 
onset of wetting occuring over such a narrow range of T around T c . 

To be explicit we consider the SU(4) calculation with L = 4. As we have seen in Fig. |3] 
the actual value of drops well below the one-loop value in eqn (|M)l once T < 2T C and it 
collapses to a value consistent with zero as T — > T c . Nonetheless we have seen that neither 
of the A;- walls shows any sign of splitting into two cd walls even at T ~ 1.02T C . This is also 
true even closer to T c , as shown in Table but the relevant runs are much shorter there. 
Consider first the k = 2 wall. Consider a fluctuation just above T c in which it breaks up 
into two cd walls separated by a region of confining phase. For this region of confining phase 
to be clearly distinguishable, it needs to be of some minimal size, presumably 0(1/T C ) long. 
Assuming (as we have argued above) that cr^T^ 2 ~ 2a c d at and near T c , the probability of 
such a fluctuation will simply be due to the free energy cost, AF^, of the region of confining 
phase. Approximating this by the latent heat, L^, multiplied by the area A and the length 
1/T C , one has a suppression 
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Plugging in the value for L h /T^ in eqn(j3SJ) and L x = L 2 = 5L we see that the range of T/T c 
in which such wetting fluctuations have a significant probability is very small and on the order 
of 

T — T 1 

- < — . (38) 

T c ~ 140 V ; 

Consider now a k — 1 wall in our SU(4) calculation. We do not know enough about its 
dynamics to analyse the situation precisely, but we can outline a plausible scenario within 
which we can see why the splitting takes place very close to T c . (The splitting essentially 
refers to the tunnelling exemplified in Fig ^] and FiglT4l) If the k = 2 wall splits into two cd 
walls at T = T c we have a^ N ^ 2 (T c ) = 2a c d(T c ) and and if we also have Casimir Scaling then 
<t^ 1 (T c ) = 4/3 x a c d(T c ). and the k = 1 wall will not split at T = T c but only at a lower 
T where the free energy advantage of a layer of confining phase provides the required extra 
2/3 x <J cd (T c ) needed for splitting. Assuming a region of confining phase of length 0{1/T C ) 
and neglecting the variation with T of both cr^ 1 and a c d leads to an estimate for the T where 
the k = 1 wall splits of 

^---^-0.013 (39) 

using eqn(|HH|) and eqn (j33|) . This is a very crude estimate, of course, but it serves to illustrate 
how the large ratio of the energy scales associated with L h /T c and a cd mean that the difference 
between the values of T at which the k = 2 and k = 1 walls break up is, assuming Casimir 
Scaling, very small and on the order of the resolution of our calculations. That is to say, our 
calculations do not in fact provide evidence against Casimir Scaling persisting at all T (until 
the walls break up). A similar argument will explain why the process of breaking up occurs 
over a very narrow range of T around T c . It also implies that any attempt to verify dynamical 
'perfect wetting' must have an extremely fine resolution in T near T c , and particularly so for 
SU(3) where the ratio Lh/T c a c d is even larger. 

We have ignored the fact that both a c d and a^ x will vary with T. This variation, as T 
moves away from T c , is likely to be driven by the growing free energy difference between the 
two phases and so is likely again to be on the order of L h /T c x (T — T c )/T c , in which case we 
will obtain a result similar to the one above. Of course, as we move away from T c either way, 
we will eventually encounter the spinodal point at which one or both types of wall ceases to 
exist. 

We briefly return to the question why a cd {T c ) is so small compared to L h /T c . At T = T c 
there is a cancellation of differences in entropy and internal energies between the confined and 
deconfined phases. Nonetheless one might have expected the 'hills' and 'passes' in the effective 
potential atT = T c to have a natural scale of O(Lh)- This would imply a value for <r^(T c ) and 
Vcd{T c ) that is 0(Lh/T c ), which is consistent with low-order perturbation theory, but not with 
what we actually find to be the case. Clearly our expectation is mistaken: the 'passes' are 
also strongly suppressed. As we move a little away from T c the free energy density difference 
between the phases increases roughly like ~ Lh\l — T/T c \. The behaviour of the domain walls 
suggests that the whole effective potential, not just the wrong vacuum, scales something like 
this: a cd (T c ) and <r^(T c ) are small because the height of the 'pass' which the tunnelling path 



21 



traverses is also small - presumably the same dynamics as the free energy cancellation. As we 
move away from T c the height of the pass also grows roughly as ~ L^l — T/T c \. This suggests 
that crj^(T) will grow in this way as T increases above T c - qualitatively in accord with what 
we see in FigOl 

This might also suggest a scenario in which both a c d(T) and cr^(T) have a minimum at 
T = T c but this requires the effective potential to be very similar in shape above and below 
T c (with just the depth of the minima reversed) which hardly seems likely. 

The fact that the fc-walls do break up only at T ~ T c does tell us that cr^ cannot be very 
much larger than 2a c d once we are very near to T = T c . This confirms the strong suppression 
of Cjy indicated in Fig. [TJ We note that this appears to be an effect that survives the N —>■ oo 
limit and provides further evidence that even in that limit the deconfined phase is far from 
being a weakly coupled gluon plasma in the region from T c to several times T c . This fits in 
with the observation that the SU(3) pressure anomaly is also a large- N effect |2*7j . 

3.4 't Hooft string condensation 

The fact that cr^(T) decreases rapidly towards a value that is close to zero, as we decrease T 
towards T c , tells us that the spatial 't Hooft string tension, a fc (T) = <r^ / (T)/T, also decreases. 

Since the deconfining transition is robustly first order (for larger N) we can use its metasta- 
bility to discuss what happens in the deconfined phase for at least some range of temperatures 
T < T c . Doing so, the observed rapid decrease suggests that &k = at some temperature 
T H that is close to and below T c . This suggests that as we decrease T in the deconfined 
phase, then at T = T H we would encounter a second-order phase transition driven by the 
condensation of spatial 't Hooft strings - if it were not for the intervening presence of the first 
order transition to the (more-or-less) separate confining ensemble of energy eigenstates of the 
theory. 

This is reminiscent of the expectation that as one increases T within the confining phase, 
there would be a second order Hagedorn-like string condensation transition at T = Th, if it 
were not, once again, for the intervening presence of the first order transition to the (more-or- 
less) separate 'gluon gas' ensemble of energy eigenstates of the theory. Indeed one can use the 
metastability of this first order transition to obtain direct evidence for this would-be second 
order transition and to estimate its location at T = Th > T c , as has been done in [2*%] . 

The condensation of confining strings as T — > T^, leads to deconfinement, but without 
reference to the 'gluon plasma'. Similarly one can see how the condensation of (spatial) 'tHooft 
strings as T — > Tj T , leads to confinement, but without reference to the usual confining phase. 
As the domain wall tension goes to zero, the vacuum will break down into a gas of bubbles 
containing different Zn phases so that the Polyakov loop, l p , is zero. 

All this hints at some kind of a duality between the confined and deconfined phases. 
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4 Conclusions 



We have shown in this paper that the surface tensions of the domain walls separating different 
deconfined phases in the Euclidean formulation of finite-T SU(iV) satisfy Casimir Scaling, i.e. 

a k w = k(N - k) f(g 2 (T)N, T/T c ) T 3 , (40) 

to a good approximation. We have shown this to be the case not only at very high T where 
our calculations merely reproduce the perturbative predictions PUE], but all the way down 
to T ~ 1.02T C where, not surprisingly, both the wall profile and the magnitude of the surface 
tension are very far from perturbation theory. It appears that the Casimir Scaling predicted 
by low-order perturbation theory is in fact a much more general property of these domain 
walls and, equivalently, the spatial 't Hooft string tension. 

Our numerical results are for da^y/dT rather than for cr^ itself, but we can interpolate 
the former and integrate it so as to estimate the T dependence of the latter. We have too 
few values of da^y / dT to be able to do so with any precision, however the striking qualitative 
feature is that as we reduce T below, say, 2T C the surface tension becomes strongly suppressed 
relative to its perturbative value, reaching values close to zero at T ~ T c . 

We have also performed calculations near and passing through T c . Here, in a very narrow 
range of temperatures around T c , we see all the k- walls splitting into confined-deconfined 
walls. The fact that the Polyakov loop profiles of k = N/2 domain walls are found to pass 
through zero, at any T, (almost) inevitably leads to the 'perfect wetting' at T = T c of such 
walls. The reason that all the fc-walls 'wet' at values of T very close to T c , is intimately linked 
with the fact that the latent heat, which provides a measure of how the confining-deconfining 
free energy difference increase with T — T c , is very much larger than the surface tension, a c d, 
of the confined-deconfined domain wall. This means that to test ideas of 'perfect wetting' for 
kj^N/2 domain walls, e.g. the usual domain walls in SU(3), will require calculations with a 
very fine resolution in T - certainly finer than ours, and quite possibly finer than any previous 
calculations. 

There are a number of ways in which our calculations could be usefully improved upon. 
Firstly, the comparison with perturbation theory would be under much better control if we 
were to have two-loop calculations of cx(y in the various coupling schemes that we have tried. 
Turning to our lattice calculations, while these do allow us some control of the continuum limit 
and the large- N limit, it is clear that in both respects our calculations could be significantly, 
and usefully, improved upon. Calculations at more values of T would transform our very crude 
reconstruction of the value of cr^(T) into something much more controlled and credible. In 
addition it would be very useful to analyse the Polyakov loop profiles, and hence the effective 
action, not just in terms of the trace, but in terms of the fuller range of gauge invariant 
variables provided by the eigenvalues of the Polyakov loop. This would allow us to transform 
our plausible inferences into strict derivations (or the opposite, as the case may be). 

We have noted that the rapid decrease of both o"{y(T) and the spatial 't Hooft string 
tension, &k(T) = cr(^(T)/T, as we decrease T towards T c , suggests the presence of what would 
be a second order phase transition a little below T c if it were not for the intervening presence 
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of the first order transition. This is reminiscent of the conjecture that as one increases T 
within the confining phase, the condensation of confining strings would occur, a little above 
T c , if it were not, once again, for the intervening presence of the first order transition. Just as 
has been done in the latter case [2H] it would be interesting to use the met ast ability of the first 
order transition to go beyond T c for direct evidence of the second order transition. It would 
also be very interesting to perform the sort of calculations we have done, but more densely 
in T, so as to determine the functional form of <7fc(T) much more reliably and precisely. This 
would help us to explore the possibility of a duality between the behaviours as T — > Tjj and 
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Appendix A 

When T T c the coupling g 2 {T) becomes large and our choice of the mean-field improved 
coupling evaluated at a = 1/LqT, i.e. gj{a), may become inadequate. To estimate the range 
of uncertainty introduced by this choice, we calculate the perturbative expression for the wall 
tension in Section EH with a 'good' coupling at the scale T _1 = aL . 

We need a coupling that has been explicitly calculated to length scales that are quite large. 
For this purpose we shall use the calculation in [TSj. This calculation is for SU(3) and we shall 
therefore perform our analysis for iV = 3. The coupling in ^J] has been extrapolated to the 
continuum limit, but we shall apply it to a finite a, assuming the corrections are insignificant 
for our purposes. 

The running of the coupling g in is given by 

= m = ~ bof ~ hl ~ 9 " ~ ^~ gl (41) 

where 

6 = H(4tt)- 2 , h = 102(4tt)- 4 , b e 2 ff = 1.5(8) (47r)~ 3 . (42) 

(Here is a fitted effective coefficient that subsumes higher order corrections and is not 
intended to be an estimate of the coefficient 62 of the /^-function in this coupling scheme.) The 
normalisation of the coupling is then fixed by its value at the largest length scale at which it 
is calculated 

f(l m )=3A8 ; l m = Aa(/3 = 5.904). (43) 

This length scale is l m ~ 2/ 3T C which is not quite as large as we would like but is large enough 
that extrapolating it to I = 1/T C might be qualitatively reliable. 
Defining the reduced length 

1 = 1 (44) 
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we solve eqns ()41l42|) numerically, with the boundary condition of eqn ()43|) . so obtaining g 2 (j) 
for any desired value of I. The quantity that appears in eqn (jl6|) and which we now need to 
evaluate is 

l_dm_d[ 

d/3g(i) g*(t) dl df3' 1 ' 

We will perform the calculation for a lattice with L = 4 and with varying T = l/a(/3)L . So 
we evaluate the coupling at 

l m L a((3 ) ay/a((3 ) ' 

We now use the expression for a^fa{(3) as a function of (3 that is provided in Table 5 of jHj 
for the range 5.6925 < (3 < 6.338, to evaluate both I in eqn (|47)|) and dl/d/3 in eqn (|43|) in 
this range of j3. (Like any multi-parameter interpolation, this cannot be assumed to be useful 
outside the range in which it was fitted.) 

We now have all the ingredients to calculate each of the factors on the rhs of eqn ()45|) . The 
first factor comes by plugging in the calculated value of / into the numerically determined 
g 2 (l). The second factor is just the /3-function in eqn(|4T|) evaluated with this g 2 (l). To obtain 
the third factor one differentiates the interpolating formula in Table 5 of [3] and evaluates it 
at the desired value of /3, and normalises it to the formula at (3q = 5.904. Since T c = l/4a(/3 c ) 
for (5 C = 5.6924(2) [2H] and this lies within our interpolation range, we are able to perform this 
calculation down to T = T c . To transform this coupling to one for SU(4) we simply assume 
that g 2 N is fixed at fixed T. 

We use the above in the formulae of Section 12.41 and show in FigQ the resulting value for 
ASw I 'L\L2 of the 2-loop perturbative wall tension using the coupling g 2 (T~ l = aLo). We also 
show the values obtained with the mean-field improved coupling, g](a), and the unimproved 
bare coupling, g 2 (a). Finally we show the values obtained from the lattice simulation. We 
observe that the more realistic the coupling, the closer the perturbative value moves towards 
the lattice values. However a substantial discrepancy remains, strongly suggesting at T ~ T c 
the domain wall can no longer be usefully described within perturbation theory. 
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Table 1: 
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Table 2: 
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Table 3: 
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Table 4: 



N 


aT 


k 


Prediction 

X _L V Ally. U ± V_y ± J_ 


AS* IL„L, 

w 1 ±J x ±J y 


AS k /AS 1 


CS 


TIT, 


2 


25 


1 


0.1042 


3(189(841 






1 88 


"i 

O 


95 


i 

± 


10^9 


1 ^80(40) 






1 09 




25 


1 


0936 


271 6(67) 






1 88 


A 
*± 


n 9^ 
u.zo 


1 
1 


U.UoU ( 


1 ri ^f^n^ 






1 09 






9 


n 11 ^7 

U. J. ±0 t 


9 1 41 f^n) 


1 ^9fi(9Q) 




1 D9 


4 


0.25 


i 


0.0791 


2200(56) 






1.88 






2 


0.1055 


0.2957(57) 


1.344(31) 


1.333 


1.88 


4 


0.25 


1 


0.0421 


0.0423(16) 






~ 1000 






2 


0.0561 


0.0567(14) 


1.341(49) 


1.333 


~ 1000 


4 


0.20 


1 


0.0514 


0.739(27) 






1.02 






2 


0.0685 


0.982(28) 


1.329(56) 


1.333 


1.02 


4 


0.20 


1 


0.0467 


0.1292(54) 






1.88 






2 


0.0622 


0.1743(60) 


1.350(57) 


1.333 


1.88 


6 


0.20 


1 


0.0381 


0.488(25) 






1.02 






2 


0.0610 


0.847(27) 


1.74(9) 


1.60 


1.02 






3 


0.0687 


1.015(29) 


2.08(10) 


1.80 


1.02 


6 


0.20 


1 


0.0345 


0.098(9) 






1.88 






2 


0.0552 


0.159(10) 


1.62(14) 


1.60 


1.88 






3 


0.0621 


0.171(10) 


1.74(15) 


1.80 


1.88 



Table 5: 



SU(4) 


P 


T/T c 


k = 1 


k = 2 


10.497 


1.016 


no 


no 


10.494 


1.012 


no 


no 


10.491 


1.007 


mixed 


flattens 


10.488 


1.002 


mixed 


mixed 


10.485 


0.998 


splits 


splits 


10.482 


0.993 


splits 


splits 



Table 6: 
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SU(6) 


P 


T/T c 


k = 1 


k = 2 


k = 3 


24.532 


1.013 


no 


no 


no 


24.526 


1.009 


no 


no 


no 


24.520 


1.004 


no 


via origin 


flattens 


24.514 


1.000 


no 


via origin 


flattens 


24.508 


0.996 


splits 


splits 


splits 


24.502 


0.991 


splits 


splits 


splits 



Table 7: 




Figure 1: Action per unit area of the k = 1 domain wall in SU(4) with aT = 0.25. Monte 
Carlo values, +, compared with perturbation theory based on verious couplings: g 2 (a), solid 
line, g](a), long dashed line, g| F (T), short dashed line. The dotted line is the interpolation 
in eqn(J3!U). 
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Figure 2: The running coupling g 2 {T) obtained from the lattice bare coupling (solid line), 
the mean field improved bare coupling (dashed line), and the Schrodinger functional coupling 
(short dashed line). All for SU(4) and aT = 0.25. 
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Figure 3: Surface tension in units of T, using the interpolation shown in Fig. 21 For comparison 
we show the 2-loop perturbative result using the mean field improved coupling. All for the 
k = 1 wall in SU(4). 
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Figure 4: The values in the complex plane taken by the Polyakov loop for 23 e [0, L — 1], for 
the k — 1 wall in SU(4) with aT = 0.25 and T ~ 1.88T C . Dashed line is one-loop perturbation 
theory. 
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Figure 5: The values in the complex plane taken by the Polyakov loop for x 3 G [0, L — 1], for 
the k — 1 wall in SU(4) with aT = 0.25 and T ~ 1.02T C . Dashed line is one-loop perturbation 
theory. 
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Figure 6: The values in the complex plane taken by the Polyakov loop for x% e [0, L — 1] for 
a sample of individual k = 2 walls in SU(4) with aT = 0.25 and T ~ 1.88T C . 
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Figure 7: The values in the complex plane taken by the Polyakov loop for x 3 G [0, L — 1] for 
a sample of individual k = 2 walls in SU(4) with aT = 0.25 and T ~ 1.02T C . 
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Figure 8: The values in the complex plane taken by the Polyakov loop for x 3 G [0, L — 1] for 
a sample of individual k — 1 walls in SU(4) with aT = 0.25 and T ~ 1.88T C . Solid line is the 
average. 
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Figure 9: The values in the complex plane taken by the Polyakov loop for £3 e [0, L — 1] for 
a sample of individual k — 1 walls in SU(4) with aT = 0.25 and T ~ 1.02T C . Solid line is the 
average. 
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Figure 10: Solid line: the average of the real part of the Polyakov loop as a function of z 
for the Lq = 4 SU(4) calculation at T ~ 1.88T C . One loop perturbative predictions using for 
g 2 {T) the lattice bare coupling (long dash), the mean field improved lattice coupling (short 
dash), the Schrodinger functional coupling (dot) and just a best fit (dash-dot). 
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Figure 11: The average of the modulus of the Polyakov loop as a function of z for the 
L = 4 SU(4) calculation at T ~ 1.02T C . One loop perturbative predictions using for g 2 (T) 
the lattice bare coupling (long dash), the mean field improved lattice coupling (short dash), 
the Schrodinger functional coupling (dot) and just a best fit (dash-dot). 
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Figure 12: The average of the modulus of the Polyakov loop as a function of z for the 
L = 4 SU(4) calculation at T ~ 1.88T C . One loop perturbative predictions using for g 2 {T) 
the lattice bare coupling (long dash), the mean field improved lattice coupling (short dash), 
the Schrodinger functional coupling (dot) and just a best fit (dash-dot). 
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Figure 13: The values in the complex plane taken by the Polyakov loop for rr 3 G [0, L — 1] for 
a sample of k — 1 walls in SU(4) at (3 C = 10.491, i.e T ~ 1.006T C . The lines represent averages 
over subsequences of 100 sweeps, taken over the first 1000 sweeps. 
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Figure 14: As in the previous figure, but for a later sequence of 1000 sweeps. 



38 



